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This movie illustrates the recent numerical and exper- 
imental discovery of a new state of turbulence in dilute 
polymer solutions, Elasto-Inertial Turbulence (hereafter 
referred to as EIT) reported in [HQ- EIT is char- 
acterized by a chaotic flow state in which both incr- 
tial and elastic effects, whose relative contributions vary 
with the Reynolds number, control the flow dynamics. 
EIT explains the phenomenon of early turbulence Q, 
which describes the onset of turbulence in the presence 
of diluted polymer additives at Reynolds numbers signif- 
icantly smaller than in the absence of polymers. EIT also 
offers a new perspective on the asymptotic drag reduction 
state in large Reynolds number polymer flows, maximum 
drag reduction or MDR Q in the limit of very large elas- 
ticity of the polymer flow. Indeed polymer additives are 
known for producing upward of 80% of drag reduction in 
turbulent wall-bounded flows through a strong alteration 
and reduction of the turbulent activity Q ■ 

In the present movie, a periodic channel flow with care- 
fully chosen initial conditions demonstrate the existence 
of EIT over a large range of Reynolds numbers. As shown 
in [l[ , the simulated evolution of the friction factor as a 
function of the Reynolds numbers matches a comparable 
experiment in a pipe flow. 

The drag reducing mechanism is caused by an increase 
of the (extcnsional) viscosity in cxtensional upwash and 
downwash flows generated by quasi-streamwise vortices 
[H, 0] , thereby creating a negative torque on these near- 
wall vortices Q . @ demonstrated that polymers re-inject 
part of the energy accumulated in high speed streaks, 
regions of locally high speed flow in the near-wall re- 
gion elongated in the direction of the flow. In inertia-less 
flows with curved streamlines, Q demonstrated the ex- 
istence of strong non-linear mixing supported by clastic 
turbulence, a state of saturated dynamical interactions 
between stretched polymer molecules and the base flow 
that causes the stretching. EIT exists by either creating 
its own extensional flow patterns, as the movie demon- 
strates in subcritical channel flows and high elasticity 
flows. 

Channel flow simulations are performed in a cartesian 
domain, where x, y and z are the streamwise, wall-normal 
and spanwisc directions, respectively. For a polymer so- 



lution, the flow transport equations are the conservation 
of mass, V • u = 0, where u is the velocity vector, and 
transport of momentum: 

<9 t u+(u- V)u = -Vp+-^V 2 u+i— -V-T. (1) 

Re Re 

The Reynolds number is based on the bulk velocity Ub 
and the full channel height H = 2h, Re = UbH/v. 
The parameter ft is the ratio of solvent viscosity to 
the zero-shear viscosity of the polymer solution and af- 
fects both the viscous stress and polymer stress terms in 
Eq. |T]). The polymer stress tensor T is computed us- 
ing the FENE-P (Finite Elastic Non-linear Extcnsibility- 
Peterlin) model 0: 
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where the tensor C is the local conformation tensor of the 
polymer solution and I is the unit tensor. The properties 
of the polymer solution are /3, the maximum polymer ex- 
tension L, and the Wcissenberg number Wi based on the 
solution relaxation time A and the flow time scale rele- 
vant to the dynamics of interest. Here Wi is based on the 
wall shear-rate 7 of the initial laminar flow at each Re, 
hence Wi = A7. The FENE-P model assumes that poly- 
mers may be represented by a pair of beads connected 
by a nonlinear spring defined by the end-to-end vector 
q. The conformation tensor is the phase-average of the 
tensorial product of the end-to-end vector q with itself, 
C = (q ® q) , whose transport equation is 

d t C + (u- V)C = C(Vu) + (Vu) T C-T . (3) 

On the right hand side of Eq. ([3]), the first two terms are 
responsible for the stretching of polymers by hydrody- 
namic forces, whereas the third term models the internal 
energy that tends to bring stretched polymers to their 
least energetic state (coiled). 

Eqs. ([T][3]) are solved using finite differences on a stag- 
gered grid and a semi-implicit time advancement scheme 
described elsewhere 111. A scries of simulations was car- 



ried out for Reynolds numbers ranging from 1000 to 6000. 
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A thorough resolution study led to choose a domain size 
of 10H x H x 5H with 256 x 151 x 256 computational 
nodes. All results discussed here have been verified on 
domains with a factor 2 in horizontal dimensions and res- 
olution in each directions. The CFL number was set to 
0.15 to guarantee the boundedness of C. 

The protocol for our simulations was designed to mimic 
the perturbed experimental setup of [l[ within the lim- 
itation inherent to the DNS boundary conditions. For 
any flow, Newtonian or polymeric, the initial flow and 
polymer fields are first equilibrated to the laminar state 
corresponding to the desired Re. A perturbation is then 
introduced over a short duration, in the form of blow- 
ing and suction velocity on both walls, over which white 
noise of prescribed intensity is introduced. The velocity 
pattern is periodic in x and z: 



,(x,z,t) = H{t) 
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where A is the amplitude, L x and L z are the horizontal 
domain dimensions, and e(t) is the random noise. The 
total duration of the perturbation is 0.5/i/t/f,, of which 
the first and last 10% correspond to a gradual increase / 
decrease through a smooth step function H(t). Choosing 
A = Q.09U b and the RMS of e at 0.005C4 causes the 



Newtonian flow to transition at Re = 6000. 

All polymeric simulations are performed with L = 200 
and (3 = 0.9. Two Wcissenberg numbers are considered, 
Wi = 100 and 700. The former is consistent with previ- 
ous simulations of MDR @, d [H. [H's theory, based 
on infinite Re and Wi, motivates the second, as an ex- 
ploration of the effects of very large elasticity, or Wi, on 
the flow. 

This work was performed during the 2012 Center for 
Turbulence Research Summer Program and will be pub- 
lished in 1151. 
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